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ASYMPTOTIC BEHAVIOR FOR A NONLOCAL DIFFUSION EQUATION 
IN EXTERIOR DOMAINS: THE CRITICAL TWO-DIMENSIONAL CASE 

C. CORTAZAR, M. ELGUETA, F. QUIROS AND N. WOLANSKI 


Abstract. We study the long time behavior of bounded, integrable solutions to a nonlocal 
diffusion equation, dtu = J *u — u, where J is a smooth, radially symmetric kernel with support 
Bd{ 0) C R 2 . The problem is set in an exterior two-dimensional domain which excludes a hole 
H, and with zero Dirichlet data on H. In the far field scale, £1 < |a:|t _1 ^ 2 < £2 with £ 1,^2 > 0, 
the scaled function \ogtu{x,t) behaves as a multiple of the fundamental solution for the local 
heat equation with a certain diffusivity determined by J. The proportionality constant, which 
characterizes the first non-trivial term in the asymptotic behavior of the mass, is given by means 
of the asymptotic ‘logarithmic momentum’ of the solution, limt_»oo f R2 u(x, t) log \x\dx. This 
asymptotic quantity can be easily computed in terms of the initial data. In the near field scale, 
|*| < W 2 /i(i) with limt-joo h(t) = 0, the scaled function t(logt) 2 M(x,t)/log |x| converges to a 
multiple of 4>(x )/ log |a;|, where <j) is the unique stationary solution of the problem that behaves as 
log |a;| when \x\ —► 00 . The proportionality constant is obtained through a matching procedure 
with the far field limit. Finally, in the very far field, |a;| > W 2 <;(t) with g(t) —> 00 , the solution 
is proved to be of order o((t log t)” 1 ). 


1. Introduction 


Let B C M. N be a non-empty bounded open set, which is assumed, without loss of generality, 
to satisfy 

( H h ) B 2 (0) CBC Bn(0), n G (2, 00 ). 


We do not assume B to be connected, so it may represent one or several holes in an otherwise 
homogeneous medium. Our goal is to describe the long time behavior of solutions to the nonlocal 
diffusion problem 


(P) 


dtu(x , t) = Lu(x, t) in (K N \H)xl 

< u(x,t ) = 0 in B x M + , 

u(x, 0) = uq(x) in M n , 


in the critical two-dimensional case N = 2, thus completing the study for other dimensions 
performed by the authors in [7J [8, [9]. The nonlocal operator L is defined by Lg := J * g — g, 
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where the kernel J is assumed to satisfy 

(. Hj ) J € C 2 (M iV ) radially symmetric, J > 0 if \x\ < d, J = 0 if |x| > d, / J = 1. 

Jr n 

Diffusion models of this kind have been widely used to model the dispersal of a species taking 
into account long-range effects [31 El CEO], and also to describe phase transitions [DEIS] and 
image enhancement na¬ 
if the initial data uq are integrable and identically zero in the hole H, problem m has 
a unique solution u € C([0, oo); L 1 (M Ar )). This is easily proved using Banach’s fixed point 
theorem; see [7j- However, in order to prove our asymptotic results for N = 2, we will need 
further hypotheses on uq. In the sequel we assume 

(Ho) uq > 0, uq = 0 in H, uq € L°°(R 2 ) n T 1 (M 2 , log \x\ dx), 


plus some extra control of the growth of uq at infinity, 

(Hi) uq G L 1 (M 2 , |x| 2 dx). 

An easy modification of the existence proof shows that if QHqD 


(1.1) u > 0, 


/ 

J R s 


u(x, t ) log |x| dx < oo 


for every t > 0, 


Moreover, if u$ satisfies in addition flHiD , then u(-,t ) € L 


holds, then 

II w IIl°°(r 2 xr + ) < II' u oIIl°°(r 2 )- 

x\ 2 dx) for every t > 0. 


Comparison with the case without holes. In the absence of holes, H = 0, the mass 
M(t) = f RN u(x, t ) dx is conserved. If the initial data are not only integrable, but also bounded, 
the solution to m behaves for large times as the solution, v, to the local heat equation 

(1.2) dtv = qAv in x K + , q := — [ \z\ 2 J(z)dz, 

2 N J r n 

with initial u(-,0) = uq ; see [6][13i. More precisely, 

lim t N ^ 2 max I u(x,t) — v(x,t) I = 0. 

oo 

Hence, the asymptotic behavior of u can be described in terms of the fundamental (self-similar) 
solution of m, 

(1.3) r q (x,i) = t~ N/2 U q , U q (y) = (47rq) _JV / 2 e _J ^~ . 

Indeed, in self-similar variables we have convergence towards the stationary profile MU q , where 

= f R JV U 0 , 

lim max \t N ^ 2 u(yt 1 ^ 2 , t) — MUJy)\ = 0. 

OO y £RiV 

Thus, there is an asymptotic symmetrization: no matter whether the initial datum is radial or 
not, the large time behavior of u is given by a radial profile, which, of course, has the same mass 
as the datum. 

The presence of holes introduces a technical difficulty, since Fourier transforms, which were 
the main tool in mm, can no longer be used. But there are differences of a more fundamental 
nature. On the one hand, mass is not conserved. On the other hand, the presence of the hole 
breaks (in general) the symmetry of the spatial domain, and an asymptotic symmetrization is 
no longer possible. It turns out that the asymptotic behavior depends strongly on the spatial 
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dimension. It is already known that there is a big difference between the case of high dimensions, 
N > 3, considered in [7], and the one-dimensional case studied in M- This paper is devoted to 
the intermediate, critical two-dimensional case. An analogous study for the local heat equation 
in dimensions N > 2 was performed in j T2| . 

Stationary SOLUTIONS. A main difference between the various dimensions has to do with 
stationary solutions of the problem, that is, functions cj) satisfying 

(1.4) L(j) = 0 in R n \ H, <j> = 0 in U. 

When N > 3, there is a unique solution of this kind approaching the constant 1 at infinity. 
Such a solution does not exist for low dimensions, N = 1,2. Nevertheless, for N = 1 there are 
stationary solutions that behave linearly at infinity. More precisely, given constants b± > 0, 
there is a unique solution to (11.41) satisfying 

(1.5) {4>( x ) ~ max{6 + x, —b~x}) € L°°(1R). 

In the critical two dimensional case there is a stationary solution with a logarithmic behavior 
at infinity, 

(1.6) (</>(x) - log |x|) € L°°(M 2 \ H). 

The construction of such solution, which is unique, is quite involved, and will be the subject of 
Section [2J 


Conservation laws and non-trivial asymptotic quantities. Stationary solutions 
play an important role in our analysis. Indeed, though mass is not conserved, under adequate 
conditions on the initial data, there is in all dimensions a conservation law of the form 

/ u(x, t)(J)(x) dx = constant 
Jr n 

for functions <j) satisfying ED and having for each dimension the right behavior at infinity 
specified above. In the particular case N = 2 we need uq £ L 1 (R 2 ,log |x| dx). 

If the initial data are also bounded, this conservation law yields on the one hand the large 
time behavior for the mass, 


M{t) -A- Ml := f > 0 if A r > 3, 

Jrn 

( Oit- 1 / 2 ) if N = 1, 

M(t) = { 

^((logf)- 1 ) if IV = 2, 
and on the other hand a global size estimate, 


u {'i t) lloo 


0(t~ N / 2 ), 

N > 3, 

or 1 ) 

N = 1, 

0((flogf) _1 ) 

N = 2; 


see Section [3] for the statements concerning the two-dimensional case. 
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Though in low dimensions there is not a residual asymptotic mass, the conservation laws allow 
to obtain non-trivial asymptotic quantities, which enter in the characterization of the large time 
behavior of the solutions. Thus, for IV = 1 we have a non-trivial limit for the right and left first 
momenta M±(t) = f R± u(x , t)|x| dx . Indeed, 

(1.7) Mf z (t) —> M*± := / uqc^ as t —>• oo, 

Jr 

where (j)^ are the solutions to (11.41) satisfying 

(<^ ± (x) — max{±x,0}) € L°°(M). 

In the two-dimensional case, the relevant quantity approaching a non-trivial constant is the 
logarithmic momentum , M\ og (t) := f R2 u(x, t) log |a:| dx. Indeed, 

M log (f) —> Ml := / uo4> as t —»• oo, 

Jr 2 

where 4> is the solution to m satisfying see Section [3j 

Outer LIMIT. The non-trivial limit quantities give an indication of the right scalings in 
order to obtain the asymptotics for the solution far away from the origin. Thus, for N > 3 the 
adequate scaling is the one that conserves mass, 

u x (x,t ) = X N ^ 2 u(X 1 ^ 2 x,\t). 

The scaled solution satisfies 

dtu x = L\u x for x € (M A \ 7~L X ), t > 0, 1~L X = {x : \ l ^ 2 x € H}, 
where L\ is the operator defined by 

L x tp(x) = \[ Jx(x-y)(<p(y)-(p(x))dy, J\(x) = X N/ ' 2 J(\ 1/2 x). 

Jr n 

If ip £ C'?°(M Ar ), an easy computation, which uses the symmetry of the kernel plus Taylor’s 
expansion, shows that L\ip converges uniformly to qAy? as A —>• oo, with q as in (11.21) . Hence, 
the asymptotic behavior is expected to be given by a multiple of the fundamental solution of 
the local heat equation with diffusivity q. Notice that the fundamental solution conserves mass. 
The proportionality constant should therefore be given by the limit mass, . Indeed, we have 

lim t N / 2 sup{|ix(a;, t) — MlTJx,t)\ : |x| > 5y/t} = 0 for all 6 > 0. 

t—too v 

Notice that the only effect of the holes in this outer limit for large dimensions, N > 3, is the 
loss of mass. 

For N = 1 the right scaling is the one that preserves the first momentum, and the asymptotic 
behavior is given in terms of the dipole solution to the (local) heat equation with diffusivity q, 

D q (x,t) = 5 x T q (x,t) = --^-T q {x,t). 

This special solution, which is self-similar, has 6', the derivative of the Dirac mass, as initial 
data, and preserves the first momentum. Thus, 

lim fsup{|u(x,t) + 2Ml + VJx, t)\ : x > 5y/t} = 0 for all 6 > 0, 

t—¥ OO ^ 
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with M* + as in (11.711 . A similar statement holds in sets of the form x < —St 1 / 2 , substituting 
M* + by . Notice that in this one-dimensional case the effect of the hole on the outer 

limit is more dramatic when compared to the case without holes: it changes both the rate of 
decay and the limit profile. Moreover, we also lose the symmetry of the asymptotic profile, even 
when the hole is small, compared to the support of the kernel, and hence does not disconnect 
the domain. 


As for the critical two-dimensional case, the fact that the logarithmic momentum has a non¬ 
trivial asymptotic limit does not show directly which is the right scaling. However, the rescaled 
mass, log tM(t), behaves for large times as twice the logarithmic momentum (it is here that 
we use the condition (Hi )); see Section 3J Hence, the rescaled mass approaches a non-trivial 
constant, namely 

log t M (t) -A- 2 Ml. 


Therefore, in outer regions the solution is expected to approach a fundamental solution of the 
local heat equation with variable mass 2M^/logt. As we will prove in Section [4l this is indeed 
the case. 


Theorem 1.1. Let N = 2. Assume that TL and J satisfy, respectively, Q Hj{ ) and ( Hj ). Let u 
be the solution to with an initial data uq satisfying QffoD Q/iip . Then, 


lim t log t sup 

t—too 


u(x,t)-2M; 


r q(x,t) 
log* 



for all S > 0, 


where T q is given by dm with q as in m, M£ = f R2 UQ<f, and f is the unique solution 
to (El and m- 


In this borderline dimension the existence of a hole modifies slightly the rate of decay in the 
outer limit, but it does not change the rescaled asymptotic profile -except for the proportionality 
constant-, which is still given by the (radially symmetric) profile of the fundamental solution of 
the heat equation. 

It is worth noticing that the mentioned scaling only gives a non-trivial profile if 0 < ei < 

| x\t 1 / 2 < £2 < OO, in the far field scale. In the very far field, |x|t 1//2 —>• 00, our result says 
that t\ogtu(x,t) —» 0, but not more. Analogous results hold for other dimensions in the very 
far field. 


Inner limit. What happens close to the holes? For large dimensions, N > 3 , solutions still 
decay as 0(t~ N / 2 ) in the inner region |x| < St 1 / 2 for some S > 0. If we scale the solutions 
accordingly, we get that the new variable w(x, t ) = t N / 2 u(x, t ) satisfies 


d t w 


T Nw 

Lw -\ -. 

2 1 


Thus, w is expected to converge to a stationary solution, solving m- To determine completely 
this solution we have to prescribe its behavior at infinity. Since there is an overlapping region 
between the inner and the outer developments, they can be matched, leading to 


t N/2 (u{x,t) - M^4>{x) T q (x,t)) -A- 0 as t — y 00 uniformly in W N . 


In particular, 




t N / 2 u(x, t ) 


uniformly in compact subsets of M N . 
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Notice that, as expected, the effect of the hole is more important when we are close to it. In 
this region, the asymptotic profile does not have radial symmetry, and sees more details of the 
hole through the function <f. 


In low dimensions the determination of the inner behavior is by far more involved. The main 
reason is that solutions do not decay at the same rate everywhere in sets of the form |x| < Dt 1 ' 2 . 
Thus, for IV = 1 the rates of decay depend on the ratio t 3 ^ 2 /\x\. More precisely, in the near field 
scale, |x| < t l / 2 h(t), lim^oo h(t) = 0, the scaled function t 3 ^ 2 u(x,t)/\x\ converges to a multiple 
of fi*(x)/\x\, where fi* is a solution to (II.4ft - (11.51) for some constants b±. The right choice for 
the involved constants is, again, obtained through a matching procedure with the outer limit. 
We obtain 


*3/2 

|x| 


^u(x,t) + 2 


4>*(x) 

X 


V q (x,t) 


-> 0 


as t —>• oo uniformly in R, 


where f>* is the solution to (11.41) (11.51) with b± = M* ± . In particular, 


t 3 / 2 u(x, t ) 


r(x) 

2q 3 / 2 y / 7T 


uniformly in compact subsets of R. 


Let us remark that in this case the matching is more complicated, since the overlapping region is 
not so wide as in the case of large dimensions. A main step is obtaining a super-solution giving 
the right decay rates in the whole inner region and up to the beginning of the far-held scale. 


In the two-dimensional case that is the subject of this paper, the rates of decay depend on 
the ratio i(logt) 2 /log |x|. Moreover, the scaled function t(logt) 2 n(a;,t)/log |x| converges in the 
near held scale, |x| < t 1,/2 h(t), Hindoo h(t) = 0, to a multiple of <p(x )/log |x|, where is the 
unique solution to <H2D and (USD- We finally obtain, after matching the inner and the outer 
developments, our main result, whose proof is completed in Section [5j 


Theorem 1.2. Under the assumptions of Theorem II. 11 the unique solution u to problem © 
satisfies 


( 1 . 8 ) 


tjlogt) 2 
log \x\ 


^ u(x,t ) — 4 M ( 


Tq(ff, t)4>(x) \ 

(log t ) 2 J 


-A 0 


as t —>• oo uniformly in R 2 . 


As a consequence, 


lim t(logt) 2 u(x,t) 

t—too 


M: 


7rq 


^ 4>(x) uniformly in compact subsets of R 2 . 


Remark. The proof can be easily extended to the case of initial data without sign restrictions. 
Indeed, if are the solutions with initial data {«o}± 5 then, by the linearity of the equation, 
u = u + — u~. Since M£ = / R2 {no}+ 1 ^ — the result for general data will follow from 

the results for u + and u~. Notice, however, that in the case of initial data with sign changes it 
may happen that M^ = 0. In this situation our result is not optimal (solutions decay faster), 
and we should look for a different scaling. 


Notice that the decay rate O ((f logt) : ) is seen not only in the far held scale, |x| ~ & 1 / 2 , 
f 0, but also in much smaller ‘parabolas’. Indeed, let h be such that h(t) —>• 0, th(t) —>• oo, 
X? ~ a i 0 < a < 1, as t —>• oo. Then, 


lim t log tu(x,t) 

t—too 


(1 - a) 


Ml 

27rq 


if lad 


th(t). 
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In particular, this result holds if h(t) = t "or h(t) = (logt) 7 with 7 > 0 (in which case a = 0). 

2 . The stationary problem 

The first aim of this section is to prove the existence of a unique solution cf to (11.41) and () 1.6 1) . 
Existence will follow from the fact that there exist ordered sub- and super-solutions with ade¬ 
quate behaviors at infinity. We will also obtain estimates for the first derivatives of (j), which 
will be used later, in Section [5l to determine the inner behavior. 

2.1. Sub- and super-solutions. We start by constructing sub- and super-solutions to the 
stationary problem behaving as log \x\ when |x| —» 00 . 

Lemma 2.1. Let N = 2, and assume that Li and J satisfy, respectively, ( | II h\ ) and ( | Hj D . There 
exist constants Ro , ko > 0 such that the function 

V-(x) = log(|x | 2 + k 0 ) - i? 0 , x € R 2 , 

satisfies 

—LV- <0 in R 2 , V- < 0 in TL. 

Proof. Since V— € C'°°(M 2 ) and the kernel J is radially symmetric, a simple calculation using 
Taylor’s expansion shows that, for q as in (E2D, 

LV-(x) - qAV.(x) = i J] f J(\x - y\)D^V-{0{x - yf, 

‘ \/ 3 \= 4 JR2 

for some ^ lying in the segment that joins x and y. We are using the standard multi-index 
notation. A direct computation yields 

4fc 0 . C 


AV-(x) = 


\D0V-(x)\ < 


|/3| = 4. 


(|x | 2 + k 0 ) 2 ’ 1 ' v, " y| (|x | 2 + k 0 ) 2 ' 

Hence, since we only have to take into account the values y such that \x — y\ < 1, 

C 


\D P V.(0\ < 


\fi\ = 4. 


(M 2 + fe 0 ) 2 ’ 

Therefore, there exists a constant K > 0, independent of ko and R, such that 

LV ~^ X) ~ (|x| 2 + fe 0 ) 2 ' 

The hrst part of our statement follows taking ko > K/ 4. 

Once ko is fixed, since Li is bounded, by choosing Rq large enough we get V- < 0 in Li. 


□ 


The construction of a super-solution, which we do next, is far more involved, and much more 
complicated than for N 2. 

Lemma 2.2. Assume the hypotheses of Lemma \ 2. 11 Given 0 < ro < d /4 such that .E? 2 ro ( 0 ) C LL, 
there are constants k > 0, 70 > 0, and D > 0 such that for every 7 > 70 there is a locally bounded 
super-solution V + of the stationary problem m satisfying 

[-LV+[x)>^, x£R 2 \H, 


( 2 . 1 ) 


V+{x) > 7 — 70 > 0, x € R 2 , 

k V+(x) = log(|x| — ro) + 7 , \x\ > D. 
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Proof. We will prove that there exist a > 0, k > 2, and a sequence 0 = ai > > • • • > Ofc+i 

such that 


v + (x) = < 


satisfies 

—Lv + {x) > in M 2 \ H 2 r 0 - 

Hence, the result will follow just taking V + = v + + 7 , with 7 > 70 := |afc + i|. 
For |x| > 2ro, a direct computation yields 


log(|x| - r 0 ), 

\x\ > D := 2r 0 + k^, 


cij, 

x € Tj := {D-j± <\x\ < D - (j - l)f} , j = 1, • • 

■ , k + 1 , 

&k +1 j 

x £ Ffc _)_3 .— B 2ro , 



A(log(|x| - r 0 )) = -- 


r 0 


| J D / 3 (log(|x| - r 0 ))| < 


C 


V 1 / 31 = 4 - 


M(M - r 0 ) 2 ’ 1 ' u/y 1 (|a;|_ ro )4 : 

Hence, Taylor’s expansion shows that there exists k € N, k > 2, large enough, such that 


( 2 . 2 ) log(|x| - r 0 ) - / J(x - y)log(|y| - r 0 ) 

Jr 2 


- r 0 ) dy > if |s| > 2 r 0 + k^- -d. 

2 (|x|-r 0 ) 4 2 . 


D 


In particular 
(2.3) 


— Lv + (x) > 


r 0 

2H 3 


if Ixl > D + d. 


The integer k may be chosen so that log(|x| — ro) > 0 for \x\ > D — d. Hence, if 0 = a\ > a 2 , 
for D < |x| < D + d we have, using (|2.2D , 

-Lv + (x) = log(|x| - r 0 ) - / J(x-y)\og(\y\-r 0 )dy - a 2 J(x-y)dy 

J{\x\ >D} J r 2 

ro|x| roD 


> log(|x| - ro) - / J(x-y)log(\y\-r 0 )dy> 
Jr 2 


> 


2(|x| — r 0 ) 4 2(D + d — r 0 ) 4 


CO 


We will prove that this estimate can be extended to the set 2ro < \x\ < D, so that 

— Lv + {x) > Co if 2ro < \x\ < D + d. 

From this inequality and (12.31) we conclude (12.11) immediately, just taking k > 0 small enough. 
To prove this claim we start by considering the annulus T 1 , and then proceed inductively towards 
the hole, choosing in each step one of the constants a,j. The positive constants 


£ = inf { [ J(x - y)dy : x eTj\, j 
1 JTi+i ' J 


— 2 , ■ • • ,k, 


will play an important role in the process. 

Let x € Ti. We already have 0 = a± > a 2 , but we are still free to choose a 2 , as long as it is 
negative. Then, if <23 < 0, since J has unit integral, 

-Lv + (x) = - J(x-y)log(\y\-r 0 )dy-a 2 J(x-y)dy-a 3 J(x-y)dy 

J{\x\>D\ JTo Jr 3 


I{\X\>D} 

> - log(-D + d - r 0 ) + |a 2 | J 2 > c 0 , 
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if 1 02 1 is large enough. Notice that 02 has to be strictly negative. 

We have already fixed ai and 02 , and have imposed that 03 < 0. If x € T 2 , and 04 < 0, 

4 

-Lv + {x) = a 2 - J(x - y)log{\y\-r 0 )dy + y^\ai\ J(x-y)dy 

■k\x\>D} Jr, 

> a 2 - log ^ - r 0 ^ + |o 3 | J 3 > c 0 , 

if |o 3 | is large enough. Since J 3 < 1, we have |a 3 | > \a 2 \. 

We now use induction. Let 3 < j < k— 1. Assume that we have already chosen the constants 
0 = ai > a 2 > ■ ■ ■ > a,j, and imposed cij+i < 0, so that — Lv + (x) > 0 for x € Let now 

X € T'j. If dj-1_2 < 0, 


i +2 

— Lv + (x) = dj + ^2 |di 
i=j~ 2 

if |ay_|_i| is large enough. Since J J+1 < 1 , we have \dj + \\ > \dj\. 

The final step is nearly identical. We have already fixed and have imposed that 

afc + i < 0. Let x € IV Then, 

fc+i 

-Lv+(x) = d k + ^2 l°i 

i=k—2 

if |afc+i| is large enough. Since J_ k+1 < 1, we have |ofc+i| > |ofc|- □ 


ITi 


J(x — y) dy > d k + \d k+1 \ J k+1 > c 0 , 


ITi 


J(x — y) dy > dj + |oj+i| J_ j+1 > c 0 , 


With the same ideas we can construct another super-solution, that will be of use in Section [5] 


Lemma 2.3. Assume the hypotheses of Lemma lff.il Let 0 < v < 1. There exist constants 
k > 0, 70 > 0, and D > 1 such that for every 7 > 70 there is a locally bounded super-solution 
w~f of the stationary problem (El satisfying 


(2.4) 


< w+(x) > 7 — 70 > 0, itl 2 , 

w t ( x ) = (log \x\) v + 7 ; M > D. 


2.2. Existence and uniqueness. The function cf solving (11.41) and m will be obtained as 
the limit when n tends to infinity of the solutions (f n to 

L4> n = 0 in B n ( 0) \ U, <f n = 0 in TL, cf n = \v- in B n+d ( 0) \ B n ( 0), 

where V— is the sub-solution constructed in Lemma 12.11 Existence and uniqueness for such 
problem is a consequence of pj, Lemma 3.1] 


Proposition 2.1. Let N = 2, and assume that LI and J satisfy, respectively, ( | Hy ) and ( Hj ). 
There exists a solution to El and El- 
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Proof. Let VL and V + be the sub- and super-solutions constructed in Lemmas 12.11 and 12.21 
Taking 7 > 70 large enough, we have \V- < V+ in M 2 . Thus, the comparison principle gives 

< (pn < D+ in B n . 

Moreover, by applying the comparison principle once again, we find that (f> n < cf n+ 1 in B n . 
In particular, there exists f> = limn-^ cj) n . It is easily seen that this monotone limit solves C3D 
and satisfies 

^ log(|x | 2 + k 0 ) - i R 0 = --V-( x ) < </>(x) < V+(x) = log(|x| -r 0 ) + C, x € M 2 , 
for some constant C > 0 , which implies (11.611 . □ 


Uniqueness follows from the fact that the unique bounded solution of (11,41) is (f> = 0. 

Proposition 2.2. Let N = 2, and assume that PL and J satisfy, respectively, Q Hu ) and ( Hj ). 
The only bounded solution to (HID is cj> = 0. 


Proof. The function <f £ = cj) — eV+ satisfies —L<f e > 0 in M 2 \ PL, and reaches its maximum at 
some finite point x, since by construction 4> e (x) —>• —00 as \x\ -A 00. A standard (for nonlocal 
operators) argument shows that if (j> e (x) > 0 we reach a contradiction. Indeed, if this happens, 
x € M 2 \ PL and (j> £ is constant in B^x). We can thus propagate the maximum to the whole 
connected component of R \ PL where x lies, which leads to a contradiction for points near the 
boundary of this component, at a distance of it less than d. Then, passing to the limit as e —>■ 0, 
we obtain <f < 0. The same argument applied to —cj) leads to > 0. □ 


Remark. This result is also true for N = 1; see [9|. However, it fails for N > 3. Indeed, in 
large dimensions there are stationary solutions that take a constant value, different from zero, 
at infinity; see j?]- 

2.3. Estimates for the derivatives. To prove them, we use that solves a problem of the 
form 

(2.5) dtu — Lu = / in R 2 x R + , u(x, 0) = uq(x), x € K 2 . 

By the variations of constants formula, solutions to (12.51) can be written in terms of the funda¬ 
mental solution F = F(x, t ) of the operator dt — L in the whole space, which can be decomposed 
as 

F(x,t ) = e~ t 6(x) + W(x,t), 

where <5 is the Dirac mass at the origin and W is a nonnegative smooth function defined via its 
Fourier transform, 

W(£,t) =e~* (V K)t -l) ; 


u(x,t) =e t uo(x) + 


W{x - y,t)u 0 (y)dy 



W(x — y,t — s)f(y, s ) dyds. 


see 0. Thus, 
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Therefore, estimates for solutions to ([2.511 . and in particular for cj), will follow if we have good esti¬ 
mates for the right hand side of the equation, /, and for the regular part, IT, of the fundamental 
solution. 

The asymptotic convergence of IT to the fundamental solution of the local heat equation with 
diffusivity q yields a first class of estimates. Indeed, for every multi-index f3 € N 2 there is a 
constant C such that 

(2.6) \D^W(x,t) - D^T q (x,t)\ <Ct~^, x€M 2 , t > 0; 

see [13]. Hence, in particular, 

(2.7) \D%W(x,t)\<Ct —5-, x <E M 2 , t > 0. 

These estimates give the right order of time decay, and will prove to be useful later, in Sections 0] 
and 0J However, they do not take into account the spatial structure of IT, and are quite poor 
for t « 0; hence they are not enough for our present goal. Instead, we will use the pointwise 
estimate 

Ct 

( 2 . 8 ) \D%W(x,t)\ < 1 + | x | 4+ | xeR 2 ,t> 0 , 


and the integral estimate 

(2.9) [ \D%W(x,t)\dx <Ct~ m/2 , x € K 2 , f > 0, 

J R 2 

both of them valid for all (3 € N 2 . These estimates were proved in m through a comparison 
argument, using that IT is a solution to 


( 2 . 10 ) 


dtW(x, t) - LW(x, t) = e~ t J{x) in M 2 x R. 
W{x, 0) = 0 in M 2 . 


Lemma 2.4. Assume the hypotheses in Provosition \2. 1\ Let 4> be the solution to m and (USD- 
There exists a constant C > 0 such that 


|V0(s)| < 

x 


in R 2 \H. 


Proof. The result will follow from an analogous estimate for ip(x) = 4>(x) — g(x), where g(x) = 
\ log(l + |x| 2 ). The function -0, which is bounded, is a solution to (12.51) with a right hand side 

f(x) = - X u {x){J * 4>){x) ~ Lg{x) : 

h{x) h(x) 


and initial data uq(x) = tf(x). 
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Let ip^(x) = tjj(Ax). By differentiating the equation for ip given by the variations of constants 
formula, we get for |x| > 1/2 and A large enough so that A x 


Vip x (x) = 


A 


1 - e" 4 


VW(Ax-y,t)ip(y)dy 


A(x,t) 


1 — e _t 


VW (Ax — y,t — s)(J * tp)(y ) dy ds 


0 JH 


1 — e 


-t 


VW (Ax — y,t — s)f 2 (y ) dy ds - AV/ 2 (Ax). 


C(x,t) 


T>(x,t) 


The first term is easily estimated using (12.911 with |/3| = 1, since ip is bounded, 

(~<\ +— 1/2 

|yl(x,t)| < —- —j- < C if we take t = A 2 , A > 1. 


1 — e 

As for the second term, the key is that, since 1~L is bounded, there exists Ao such that for 
y S TL and A > Ao there holds that |Ax — y | > ^A|x|. Hence, using that J * <j) is bounded in TL, 
and estimate (|2.8I) with |/1| = 1, we have 


\B(x,t)\<c\n\ 


A 


1 — e~ l A 5 |x| 5 


/ (t-s) 
Jo 


— s) ds = 


ca~H 2 

(1 — e _t )|x| 5 


< C if |x| > 1/2, A > A 0 , t = A 2 . 


In order to control the third term we decompose the integral in two parts and use the bound 
I/ 2 (a:) | < 1 which was obtained in the course of the proof of Lemma 12.11 We have, 


\C(x,t)\ < 


1 — e t 


0 J{\v\<m) 


\VW(Ax-y,t- s)\ 


dy 


(1 + \y\ 


212 


ds 


+ 


1 — e 


-t 



1 0 J{\y\>±\x |} 


Ci(x,t ) 

|VVL(Ax — y,t — s)| 


dy 


(1 + H 


212 


ds . 


C2(x,t) 

The hrst integral is easily controlled thanks to (I2.8|) with \f3\ = 1, 


Ci(x,t) < 


C A t 2 


j 

J \h 


dy 


1 


< C if |x| > -, A > 1, and t = A 


1 -e 1 A 5 |x| 5 J{\y\<^\ x \} (l + |y| 2 ) 2 2’ 

On the other hand, using the estimates (12.711 and ()2.8I) with \j3\ = 1, we get 


C 2 (x,t) < C 
< C 


A 


1 — e t 
A" 1 


t-i rt 

(t — s)-*/ 2 ds+ (t — s) ds 

Jt -1 

1 


dy 


( 1 + \y \ 2 ) 2 


ifW - 2 ’ X ~ 1 ’ * = A ' 
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We finally estimate V. To this aim we notice that d x% f '2 (x) = L J ■ Hence, since 


1 + x 


< 


C 


(1 + |x | 2 ) 2 ’ 


D 13 


Xi 


1 + x 


< 


C 


(l + |x | 2 ) 2 


, l/3| = 4, 


using Taylor’s expansion we get \d Xi f 2 {x)\ < , which implies 

C\ 


\V(x,t)\ < 


(1 + A 2 |x | 2 ) 2 


< C if |x| > 1 / 2 . 


In conclusion, |V^ A (x)| < C if \x\ > 7 and A > Ao- This yields, by taking |y| > Ao, x = y/\y\ 
and A = |y|, 

\y\\V^y)\ = \V^(x)\<C. 

Recalling the definition of ijj, we see that this implies that 

|V^(y)| < if \y\ > A 0 . 

\y\ 


The lemma is proved, since cj) € C 1 (M 2 \ T~L). 


□ 


Remark. As a corollary of Lemma 12.41 we have that 

( 2 . 11 ) \cf>(y) — 4>(x)\ < T^-r if x G M 2 \ "H, \x — y\ < d. 

\x\ 

Indeed, applying the Mean Value Theorem, 

... . .. Cd Cd 2 Cd , 1 1 ,1 1 

\(j}(y) — 4 >(x)\ < 7 - 7 - < 77- - < 7—7 if x € R" \ RL, x > 2d, x — y\ < d. 

Isl |*| — d |x| 

The estimate for |x| < 2d is immediate, since (j) is locally bounded. 


3. A CONSERVATION LAW AND SOME BY-PRODUCTS 


Solutions to m satisfy a conservation law that will allow us to prove that the mass has a 
logarithmic decay rate. This is later used to obtain the decay rate of solutions, a super-solution 
that gives the right rates of decay in inner regions, the asymptotic limit of the ‘logarithmic’- 
momentum, and bounds for the first and second momenta when flUiD holds. 


Proposition 3.1. Let N = 2. Assume that Li, J and uq satisfy l \H^ p , i ] // / [ ) 
tively. Let cj) be the solution to m and m, and u the solution to problem 


and QlLoD respec- 
dPl) . Then, 


(3.1) 



u(x, t)(f>(x) dx 



for every t > 0. 


The proof, based on Fubini’s Theorem, is entirely similar to the one for dimensions N=1 or 
N > 3; see [30J. Hence, we omit it. 

A first consequence of this conservation law is a bound for the mass at time t which in 
particular shows that the mass decays to 0 as t —>• 00 . 
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Proposition 3.2. If in addition to the assumptions of Proposition \3.1\ we have also uq € 
L°°(M 2 ), then there exist to > 0, Co > 0 such that 

(3.2) M(t) := f u(x,t)dx <-—— ift > to- 

J R 2 logt 


Proof. The solution, uc, of the Cauchy problem with initial condition uq is a super-solution of the 
Cauchy-Dirichlet problem (|Pl) . since uc(-,t ) > 0 for t > 0. Therefore, the asymptotic behavior 
of uc, proved in m, provides an estimate of the time decay rate of u: u(x, t ) < uc{x, t) < Ct 1 . 
Though this rate is not optimal, as will turn out later, it is enough to get the decay rate of the 
mass. Indeed, let cr > 0 such that log \x\ < a(f(x) for x ^ Pi (recall that f) > 0 on dTL). Then, 


M(t) < / u(x,t)dx + a 


\x\ 2 >5(t)t 


u(x, t)(j>(x ) 
log |x| 


dx < Ct 1 5t + 


a 


-- Ml 

\ log (St) 0 


if 5(t)t > V?, with 1Z as in ( | / ly | ) . Taking 5(t) = 1/logf we get that both terms on the right 
hand side of this estimate almost balance, and (13.211 follows. □ 


Remark. We will see later that log tM(t) has a nontrivial limit as t —>■ oo. Thus, the decay rate 
in Proposition 13.21 is optimal. 

With this estimate for the decay rate of the mass we can improve the decay rate of u(-,t). 

Proposition 3.3. Under the assumptions of Proposition 1 ,9. i?L there exists a constant C\ > 0 
such that 

(3.3) u(x,t) <-Q— ift>t 0 . 

tlogt 

Proof. We estimate u at time t by the solution of the Cauchy problem with initial data u(x, i/2). 
Thus, using the bound (12.71) with |/3| = 0 for W, and the bound for the mass (13.21) . we get 

u(x,t) < uc(x, t) = e^ l u{x, i/ 2 ) + / W(x — y, t)u(y, i/ 2 ) dy 

Jr 2 

< e_ *II w o||l°°(k 2 ) + Ct~ l M{t/ 2 ) < if t > i 0 . 

□ 


As we will see in Section [U the bound provided by (13.31) gives the optimal global decay rate 
for the solution. However, u decays faster than O ((ilogi) -1 ) in inner regions. Once we have at 
hand the global decay rate, this last assertion is proved by means of comparison with a suitable 
super-solution that we construct next. A similar super-solution will be used later in the study 
of the inner limit. With that application in mind, we keep a parameter 7 > 2 that will be set 
equal to 2 in the estimate of u. 


Lemma 3.1. Let 0 < a < q = ^ f R2 J{z)\z\ 2 dz, and let V+ be the super-solution (depending on 
the radius ro > 0) to the stationary problem constructed in Lemma 12.21 There exist b, T > 0 
such that the function 


V(x,t) 


T a (x,t)( V+(x) + b) 
log 7 1 


I \(x,t) 


e 4at 

47 t at 
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satisfies 


(3.4) (d t V - 


Proof. We have 


t r T ~\r\{ 1 ^ r aM) 0 + 0*0 + &) . 

V-LV)(x,t) >-{--!) - q - »*, 


o < \x\ 2 < 2 at, t > T. 


(dtV - LV)(x,t) = V+ } X \+ b (dtT a - LT a ) (x,t) 

log' t 

r a (x,t)(V+(x)+ b) r a (x,t)LV + (x) 


(3.5) 


- 7 - 


t log 7+1 1 


log 7 t 


1 


log 7 1 


J(x - y)(T a (y,t) -T a (x,t))(V+(y) - V + (x))dy . 


A(x,t) 

In order to estimate <9jT a — LT a , we observe that Taylor’s expansion yields 

LT a (x,t ) = qAT a (x,t) + -jy ^ f J(\x - y\)D^T a (^,t)(x - yf dy 


101=4 ’ 


for some £ lying in the segment that joins x and y. Hence, since dfT a = aAT a , 

(9 4 r a - LT a )(x,t) = (l - ^jd t T a (x,t) - [. J (\ x ~ y\) D x r a{^t)(x - yf dy. 

° ' \/ 3 \= 4 Jm2 

To deal with the main term, we compute 

An estimate for the error term follows from 

\D%TM,t) < ^ < CTa ^ t] if |/3| = 4, |x| 2 < 2 at. 

Thus, since 0 < a < q, for some to large enough we have the lower bound 


(3.6) ($r a -LT a )(M)>^Q-l) 


T a (x,t) CT a (x,t ) > 1 (q 


- 1 


r a (x,t) 


t > t 0 - 


t t 2 4 V a J t 

Let us now bound the integral term A(x,t). On the one hand, 

I V + {y) - V+ (x) | < log(D + d-r 0 )~ a k + 1 if x,y G B D+d { 0). 

On the other hand, by the Mean Value Theorem, 

\v+{y) - v + (x)\ = | W + (£) • (y - x)\ < ^ if x,y £ M 2 \ B D ( 0). 

|£|-r 0 \x\ — d — r Q |x| 

The difference between the values of T a at x and y is also estimated thanks to the Mean Value 
Theorem, using that 

|vr„(£,t)| = < C\x\T a (x,t) 


2 at 


t 


Gathering all this information, we obtain 
(3.7) 


\A(x,t)\ < for 4rp < \x\ 2 < 2at. 

t log ' t 
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If we use estimates (13.611 and (13.71) , together with (12.11) , in (13.51) , we get 

T a (x,t)(V + {x) + b) 


-c d . 




10 \a 


t log 7 1 


if b is large enough and t >T, with T large. 

Now we can obtain the time decay rates in all near field scales. 


□ 


Corollary 3.1. Let 0 < a < q. Under the assumptions of Provosition HPl there exist constants 
C,b,T > 0 such that 

tl(x , t) < c £^h±M±h j/N 2< 2ot , t> t. 

Proof. Once we have established the global decay rate GL 2 ), the result follows by comparison 
with a large enough multiple of the super-solution V constructed in Lemma 13.11 with 7 = 2 , 
taking T > to, since 


V(x,t ) > 


t logf 


if 2 at <\x\ 2 < 2 at + d 


for a certain constant k > 0 , and 

V(x,T)>c>0 if |x | 2 < 2 aT. 


□ 


A fourth by-product of the conservation law is the large time behavior of the “logarithmic 
momentum”. 


Corollary 3.2. Under the assumptions of Provo sition VLB. 

lim / u(x, t) log \x\ dx = Ml. 
^°° J R 2 


Proof. Using the conservation law (13.11) . together with the mass decay estimate (13.21) and the 
estimate (11.61) for the stationary solution, we immediately obtain 

r (j 

/ u(x, £)| log \x\ — 4>(x)\ < if t > to, 

J R2 


/ 

J R 5 


u(x, t ) log |x| dx — M ^ 
from where the result follows. 


< 


logf 


□ 


As mentioned in the Introduction, if uq satisfies ( jffiD , the second momentum of u(-,t) is 
finite for all times. We now obtain an estimate for its growth, which follows from the mass 
estimate m- This estimate for the second momentum will be used later to study the limit of 
the rescaled mass. 


Lemma 3.2. If, in addition to the assumptions of Provo sition HOI uo satisfies QiliD , then there 
are constants t\ > 0 and C 2 > 0 such that 

t 


M 2 (t) := / u(x, t)\x\ 2 dx < C 2 

J R 2 


log* 


(3.8) 


for t>t\. 
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Proof. There holds, 

M' 2 {t) = / dtu(x, t)\x \ 2 dx = 


Lu(x, t)\x \ 2 dx 


= / Lu(x,t)\xf dx — / / J(x — y)u(y,t)\x \ 2 dy dx 

Jr 2 Jh Jr 2 


< 


[ Lu(x, t)\x \ 2 dx = [ u(x, t)L\x \ 2 dx 
Jr 2 Jr 2 

c / u(x,t ) < 

Jm 2 


cCp 
log t 


for t> to- 


Integrating this inequality, we obtain (13.81) . for some t\ > to- 


□ 


The control of the growth of the second momentum yields in turn a control on the growth of 
the first momentum. 


Lemma 3.3. Under the assumptions of Lemma 1 3. 21 

r t 1 / 2 

(3.9) Mj (t) := / u(x,t)\x\ dx < C 3 - - fort>t\. 

Jr 2 log t 

Proof. Using (13.21) and (13.81) . we get 

Mi(t) < t 1 / 2 / u(x,t) dx +—pr / \x\ 2 u(x, t) dx < 

JlxKt 1 / 2 t ' J\x\>t'/ 2 


4. Outer limit 


Cpt 1 / 2 C 2 t 1 / 2 
log t log t 


□ 


We have now all the necessary tools to obtain the outer limit, Theorem 11. II The first step is 
to determine the asymptotic value of the rescaled mass in terms of the initial condition. 


Proposition 4.1. Under the assumptions of Theorem \1.2l 

log t / u(x, t) dx —>• 2Ml as t ^ 00 . 

Jr 2 

The result follows immediately from Corollary 13.21 combined with the following lemma, that 
states that the scaled mass behaves for large times as twice the logarithmic momentum. 


Proposition 4.2. Under the assumptions of Theorem \1.2l 

lim logt / u(x,t)dx= lim 2 / u(x,t)^og\x\ dx. 

t-Kx) J R 2 t->oo 7 m 2 

To prove this proposition we consider separately the regions {|x| 2 < tlogf} and{|x| 2 > flogt}. 
Lemma 4.1. Under the assumptions of Provosition HOI 

(log log t ' 


/ u(x,t) (lost — log Ixl 2 ) dx = O . 
{N 2 <tiogt} n J V log* 


as t —>• 00 . 
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Proof. We have 
\A(t)\ < 




u(x,t ) 


log — 


dx + 


{logT —l x l 2 —^ lo § 


u(x,t ) 


tog — 


dx. 


Notice that 


Ai(t) A2 (t) 

I |2 

log < Clogt in {4 < |x| 2 < i/logf, t > e}. Therefore, since u(x,t ) = 0 for 
\x\ < 2, using the size estimate m we get 

C log t f 


\Mt)\ < 


tlogt J{\ x \ 2 <^- t } log t 


(J 

dx = - -, t > max{to,e}. 


Using now that 


of the mass (|3.2jl . we obtain 


1 \ x \ 2 

log 


< Clog log t in {4/logt < |x| 2 < tlogt, t > e}, together with the decay 


U 2 (t)| < Cloglogt [ u(x,t)dx < ^ og °gt ^ t>max{to,e}. 
JR 2 log t 


□ 


In the case of the second region, both the integral corresponding to the mass and the one 
corresponding to the logarithmic momentum converge separately to 0. 


Lemma 4.2. Under the assumptions of Lemma \3.S\ 


logt / u(x,t)dx = 0 

J {|a:| 2 >t logt} 


1 


logt 


J {\x\ 2 >tlogt} 

Proof. First, using Lemma 13.2 
0 < log t 


u(x, t) log \x\ 2 dx = O 


1 


log tj ’ 


as t —> 00 . 


'{M 2 >il°gi} 


u(x,t)dx< —— [ u(x. t) I x 1 2 dx < ———— t>max{ti,l|. 
' ~ tlogt V - logt - 


As for the second integral, we use that the function r —> is nondecreasing for r > e. 

Hence, taking t 2 such that t 2 logt 2 = e 2 , and using again Lemma [3721 we have 

f / , \ 1 I ,2, f / .X M 2 , log(tlogt) 

/ u{x, t) log \x\ ax = I u(x,t) Ul2 ax < 

J {|rr| 2 >t logt} 


{|a;| 2 >ilogt} 


IN 2 


t logt 


■ / u(x,t) I 
Jr 2 


x\ 2 dx 


log(tlogt) t C , , 

< C 2 —:—— -—-<-—-, t > max{ti,t 2 }. 
t log t log t log t 


□ 


Once we have identified the limit of the rescaled mass, we can proceed to obtain the outer 
limit. Let us remark that, although we are dealing with the outer region, the proof requires the 
refined bound in the inner region provided by Corollary 13.11 

The idea is that logtu(-,t) ~ u(-,t) for large times in the outer region, where v = v(x,s) is 
the solution to 

(4.1) d s v — qAu = 0, x € M 2 , s > r(t) := f/(logf) 1 ^ 2 , u(cc,r(t)) = log r(t)u(x, r(t)), x € M 2 . 
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Lemma 4.3. Under the assumptions of Theorem li.il the solution u to problem [71 satisfies 
^liin tsup j|logf u(x,f) — v(x,t)\ : |x| > 6\/i| = 0 for all 5 > 0, 
where v is the solution to dm 


Proof. Since u is a solution to the Cauchy problem (12.51) with right hand side / = —Xh(J * u )> 
and initial data u(-,r(f)) at time s = r(t), then 

t{ log tu(x,t) ~ v{x,t)) = t(log t — logr(t))u(x,f) + 1 logr(i) u{x,r(t)) 

'---V-' '-----' 

A(x,t) B(x,t) 

+1 log r(t) [ (W(x-y,t-T(t))-r q (x-y,t-T(t)))u(y,T(t))dy 
Jr 2 

' -v-' 

C(x,t) 

— t log r(t) [ e~ {t ~ s) xn(x)(J *u(-,s))(x)ds 

Jr{t) 

^-v- 

T>(x,t) 

-tlogr(t) / / W(x-y,t-s)(J*u(-,s))(y)dyds. 

Jr(t) Jh 

s ’ v y 

£(x,t) 

We have, using (13.31) . the L°° bound in (11.11) , estimate m with |/31 = 0, the bound for the 
mass (13.21) . and hypotheses ( | //^ D ■ 

\A{x, t)| = |1 log log t u(x, t) < Ci loglogt/(21ogf), 

\B(x,t)\ < tlogr(i)e“( t-T ^||u 0 || i ao( R a), 


\C(x,t)\<Ctlog(r(t))(t-r{t)) 3/2 M(r(f)) < Ct(t - r(t)) 3/2 , 
T>(x,t ) = 0 if \x\ 2 > 5 2 t, t > 1Z 2 /<5 2 , 


which implies the uniform convergence to 0 of all these terms as t —>• oo. In order to check that 
also 6 converges to 0, we observe that u(z,s) < for \z\ < 1Z + d] see ( | II A ) and Corol¬ 
lary [3TTJ Therefore, ( J*u(-,s )) < in 7L Now, by using that W(x,t) < Ct/ |x| 4 , see (|2.8D . 

we get, for \x\ 2 > S 2 t with t large enough so that \x — y\ > \\x\ for every y GPL, 


\£(x,t)\ < f logr(t) 


C 


t — s 


\x\ 


r(t) Sfog 


ds < Ct fog r(t)(t - T(t)) 2 < 


C s t 


r(f)(logr(f)) 2 |x| 4 r(t) log r(f) 


0 as t —>■ oo. 


□ 


In view of Lemma 14.31 Theorem 11.11 will follow from the large time behavior for v(-,t) given 
next. Notice that the result is not immediate, since the initial data for v moves with t. 

Lemma 4.4. Under the assumptions of Theorem 11.11 the solution v to m satisfies 

lim t sup{|u(x, t) — 2 M£Tq(x, f)| : x € M 2 } = 0. 

f. —VfYl t ^ 
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Proof. We have 


t 


v(x,t) - 2M1T q (x,t) = t\v{x,t) - logT{t)M(T(t))r q (x,t - r(t))\ 

A(x,t) 

+1 log T{t)M(T(t)) \T q (x,t- r(t)) - r q (x,i)| 

^ v X 

B(x,t) 

+ t|r q (x,t)| |log T{t)M(r(t)) — 2M^| . 

'-V-' 

C(x,t) 


Notice that tr q (a;,i) € L°°( R 2 ), and 


lim £sup{|r q (x,i) — T„(x,t — r(t))| : x € M 2 } = 0. 

t—> oo 

Hence, Propositions 13.21 and 14.II imply, respectively, that B and C converge, uniformly in x, to 0 
as t -A oo. Therefore, since t — r(t) > f/2 if i is large enough, the result will follow if we prove 
that 

^lim (t — r(t)) sup { \v(x, t) — log r(t)M (r(t))r q (x, t — r(t))| : x € R 2 } =0. 

This is what we do next. 

Since v is a solution to m, its Fourier transform satisfies 

w(£,s) = v(£,T(t))e- 4n2q ( s - T ^\\ 

Moreover, D(0, r(t)) = logr(t)M(r(t)). Therefore, performing the change of variables rj = 
(s — T (*)) 1/2 «. we have, for some function R(t) that will be conveniently chosen later, 

0 - r(t)) v(x,s) - logT(t)M(r(t))r q (x,s - r(t)) 


(s — r(t)) f (w(C,r(t))-t)(0,r(t)))« 

Jr 2 

sup |t)(£,r(i)) - u(0,r(t))| [ e~ 4n2qM 
I <R(t) JR N 


< (s — 

< 


—47r 2 q-(s—-3-(i))|C| 2 


df 


i€i <m 

+2 log r(t)M(r(t)) 


dr] 
— 4 tt 2 q\v \ 2 


dr/. 


'\v\>(s-r(.t))V*R(t) 

On the other hand, using the Mean Value Theorem and (13.91) . we obtain that 

sup |D(^, r(t)) — D(0, r(t))| < 2i?(t) log r(t) f |x||n(a;, r(t))| dx < C7?(t)(r(t)) 1,/2 . 
ICI <R(t) 4r2 

Thus, taking s = t, and considering Proposition 13.21 we arrive at 

(t-T(t)) | v(x,t) - log r(t)M (r(t))r q (x, t — r(f))| 


<CR{t){T{t)) l > 2 + C [ 

J Ini 


e dr/ ^ 0 as t -A oo, 


>a —47T 2 q|77| 2 

'\V\ >(t-r(t))Vafl(t) 

if R{t)(r{t)) 1//2 —» 0 and (t — r{t)) l t 2 R(t) -A oo as t -A oo. These two conditions are fulfilled if, 
for example, 

i?(t) = f _1/2 (logf)“, 0 < a < 1/4. 

□ 
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5 . Inner limit 


The goal of this section is to complete the proof of Theorem ll.21 In view of Theorem ll.il what 
is left is to show that the limit m is valid uniformly in an inner set of the form {|x | 2 < 5t} 


for some 5 > 0. Since 


4>(x)/ log |x| 


< C in M 2 \ TL, this will follow from the next result, if we 

take into account (12.61) with |/3| = 0. 

Theorem 5.1. Under the assumptions of Theorem 1 1. 21 for every 0 < a < q we have 


(5.1) lim i(logt) 2 sup 

£—>•00 


log |x| 




ct>(x)W(x,t ) 


(logt ) 2 


x € 


\ R, |x '| 2 < 2 at > = 0 . 


The advantage of this formulation in terms of W is that it is more straightforward to apply 
the nonlocal operator L to W than to r q . 

In order to prove m we perform a comparison argument in {|x | 2 < 2at) with suitable 
barriers which are fast enough e-close to the asymptotic limit as t —>• oo, 


(5.2) 


w. 


CM) = 4. m; 

_ 


4>(x)W(x,t) 
(log t ) 2 

N/ * 

v(x,t) 


±Re(x,t), 


Re> 0, 


lim t(\ogty 

t— KX) 


sup 


Re(x,t) 
log |x| 


We start by estimating how far is v from being a solution, since this will be needed to prove 
that uif and u~ are respectively a super- and a sub-solution. 

Lemma 5.1. There exists a constant C > 0 such that 

C 

\dtv — Lv\(x,t) < -ut ifx€R 2 \T~L, Ixl 2 < 2at, t > 1. 

t z (logt) z 


Proof. A straightforward computation shows that, for x € M 2 \ TL, 

_J_ _ (dtW - LW)(x,t)</>(x) 2 cf(x)W(x,t) 

4 tV V X ' (logt ) 2 t(logt) 3 

+ n 1 .s 2 [ J(x-y)(W(y,t)-W(x,t))((j)(y)-(j)(x))dy. 

(log t) z J R 2 

' -v-' 

A(x,t) 


In order to estimate the integral term, we write it as 

A(x,t) = Vr q (x,t)-/ J(x-y){y-x)((j){y)-<t)(x))dy 

Jr 2 

+ (VW(x,t)-VT q (x,t)) ■ [ J(x-y)(y - x)(f(y) - f>(x))dy 


l J ( x -y) D x^(C,t)(y-x)^(^(y)-^(x))dy. 
2 |^ 2 7k2 


Estimating the factors involving derivatives of W by (12.61) with \fd\ = 1 and (12.71) with |/3| = 2, 
the gradient of T q by |Vr q (x, f)| < C\x\/t 2 , and the difference of the values of f> at y and x 
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by (12.111) . we get that \A(x,t)\ < C/t 2 if x £ M 2 \ T-L. Therefore, since W is a solution to 
problem (12.101) . we conclude that, for x £ M 2 \ R, \x\ 2 < 2 at, 


\d t v - Lv\(x,t ) < 4 

(. 


e t J(x)4>(x) | r) 4>(x)W(x,t) 

+ Aj — ——3 + 


c 


(logt) 2 


< C 


' t log |x| 


+ 


t(logt) 3 

1 


f 2 (logt) 2 

C 

< 


l (logt) 2 f 2 (logt) 3 t 2 (\ogt) 2 ) f 2 (logt) 2 


□ 


We now turn our attention to the correction term R £ . Notice that, in addition to having the 
required decay property specified in (15.21) . R e should be such that uf and u~ are respectively a 
super- and a sub-solution. To this aim, we need dtR e — LR e to be strictly positive, and with a 
slower decay than |<9 t v — Lv |. 

Guided by our experience with the non-critical cases N A 2, our first attempt is to take a 
function in separated variables. A good try is to pick v € (0,1), and take 


wj{x) 

C v t(\ogt) 2 ' ^ 

w v (x,t) 


c v 


sup 

xeM 2 Y h 


w+(x) 

log |x| 


where is the super-solution to the stationary problem (11.41) given by Lemma 12.31 for this 
value of v. Notice that this choice leads to the right decay for R e . But we still have to check 
whether oof is a super-solution. Hence, we have to compute dtw v — Lw v . 

Using (GUD, we get, for a certain constant k, 


(5.3) 


( d t w v - Lw u )(x,t ) = — 


w+(x) 


(1 + 2(logt) ^-f- 


f 2 (logf) 2V " ' ’ ' |x| 2 (log \x\) 2 ~ v 1 (logt) 2 ' 

Therefore, using again (EZD, we obtain that there is a constant a > 0 such that 


(d t w u - Lw u )(x,t ) > 


(/tt-cr|x| 2 (log|x|) 2 ) 2 


|x| 2 (log \x\) 2 ~ u t 2 (logf) 2 


, x € M 2 \ R, t large enough. 


Since |x| log |x| < ( St ) 1 / 2 if |x| < ^ ^ and t > nrax{(5, e}, taking 5 = k/{2g) we finally obtain, 
for some large enough T, 

Kj St 

(d t w„-L^)( x ,t) > — 2(logt)2 _„ . *€« 2 \«. N < t>r. 

and we are done in this region. Unfortunately, we can not even guarantee that w u is as super¬ 
solution in the set j — \ x \ 2 — 2at|, hence the need to look for an alternative. 

A second possibility is to take R e = KV, where V is the super-solution that was constructed, 
for a fixed 7 > 2, in Lemma 13. II and K is a positive constant to be fixed. This choice also has 
the right decay specified in (15.21) . Moreover, we see from (|3.4I) . that 


d t V -LV> 


C 

t 2 (logt)'i' _1 ’ 



St 

(log t) 2 ' 


t > T. 
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Hence, if 7 € (0,3), this alternative works in this region. However, our estimate from below for 
dtV — LV decays too fast in the set jx € K 2 \ PL : \x \ 2 < j; and is not able to control the 

contribution of v in order to make 1 uf a super-solution there. 

What will work is a combination of the two possibilities that we have considered, 


R e = R e ^K(x,t) = £w v {x,t) + KV(x,t), K > 1, e > 0. 


Indeed, since V is a super-solution, by the very same computation performed above we have 

l^j n . , ,0 St 


d t u>t - Lujt > 


x € M 2 \ PL, \x\ 2 < 


t > T. 


2St 2 (logt) 2 ~ l/ ’ ’ (logf) 2 ’ 

On the other hand, we see from (15.31) and (12.4p that, for some large enough T, 

C S 2 t 


\d t w v - Lw v \(x,t) < 


x£R 2 \R, 


2 < \ x \ ^ 2at, t >T. 


t 2 (logt) 2 ~ L '’ ~ ^ (logt ) 2 

Thus, using (13.41) . we see that there is a large enough time T > 0 such that, for every K > 1, 


O + t + \ (C\og\x\ C 

d t u e Lu e _ y t 2( logt p t 2 {\ogt ) 2 


eC \ 
t 2 {\ogt) 2 ~ u ) ’ 


x£R 2 \H: \x \ 2 < 2 at, t >T. 


Therefore, if we first choose 2 < 7 < 3, and then 0 < v < 3 — 7 , we finally obtain that uf is a 
super-solution in {x € H 2 \ R : \x\ 2 < 2 at, t > T e } for some T e > 0 large enough. Analogously, 
for every e > 0 and K > 1, co~ is a sub-solution in the same set. 


Proof of Theorem \5.1\ Let e > 0. The outer behavior, Theorem 11.11 together with (12.61) with 
|/3 1 = 0 , implies that there exists a time t £ such that 


t(logf) u(x,t) - 2 


< 2e, 


2 at < \x\ 2 < 4 at,t > t £ . 


Since -A- \ and -A- 1 uniformly in { 2 at < \x\ 2 < 4at} as t —> 00, 


u(x,t)-4M* 


(t>(x)W(x,t ) 
(log t) 2 


< 4 e 


log |x| 
f(logf ) 2 


< Ce4Ml 


4>(x)W(x,t) 
(log t) 2 


if 2 at < \x\ 2 < 4at, t > t £ . 


Thus, if e is small, 

(1 — Ce)(v(x,t) — R £j K( x ,t)) < u(x,t ) < (1 + Ce)(v(x,t) + R £) K(x,t)), 2 at < \x\ 2 < 4at, t > t £ , 
since R e = R £) k > 0. On the other hand, if we take K > 1 large enough, 

(1 - Ce){y - R £yK )( x ,i £ ) < u(x,t £ ) < (1 + C'e)(u + R £yK )(x,t £ ), x £M 2 \'H, \x \ 2 < 4at E . 

Hence, by the comparison principle we get, 

(1 — Ce){v — R £) k) < u < (1 + Ce)(v + R £ ,k) x € M 2 \ PL, \x\ 2 < 2 at, t > t £ . 
Therefore, using the decay properties of R £) k , we conclude that 


lim sup t(log t) 2 sup 

t—t OO 


log \x\ 


| u(x, t ) — u(x, t)| : x € M 2 \ R, \x\ 2 < 2 at > < Ce. 


Since e is arbitrary, we get the desired result. 


□ 
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